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We present RlCO, a code designed to compute the ionization fraction of the Universe during 
the epoch of hydrogen and helium recombination with an unprecedented combination of speed and 
accuracy. This is accomplished by training the machine learning code Pico on the calculations 
of a multi-level cosmological recombination code which self-consistently includes several physical 
processes that were neglected previously. After training, RlCO is used to fit the free electron fraction 
as a function of the cosmological parameters. While, for example at low redshifts (z < 900), 
much of the net change in the ionization fraction can be captured by lowering the hydrogen fudge 
factor in Recfast by about 3%, RlCO provides a means of effectively using the accurate ionization 
history of the full recombination code in the standard cosmological parameter estimation framework 
without the need to add new or refined fudge factors or functions to a simple recombination model. 
Within the new approach presented here it is easy to update RlCO whenever a more accurate full 
recombination code becomes available. Once trained, RlCO computes the cosmological ionization 
history with negligible fitting error in ~ 10 milliseconds, a speed-up of at least 10 6 over the full 
recombination code that was used here. Also RlCO is able to reproduce the ionization history of the 
full code to a level well below 0.1%, thereby ensuring that the theoretical power spectra of CMB 
fluctuations can be computed to sufficient accuracy and speed for analysis from upcoming CMB 
experiments like Planck. Furthermore it will enable cross-checking different recombination codes 
across cosmological parameter space, a comparison that will be very important in order to assure 
the accurate interpretation of future cosmic microwave background data. 



I. INTRODUCTION 

Planck, the third generation of satellite missions to 
study the Cosmic Microwave Background (CMB), is 
scheduled for launch later this year. It will make ac- 
curate measurements of the temperature fluctuations in 
the CMB out to I ~ 2500 as well as provide a detailed 
picture of the CMB polarization to I ~ 1500 [![. In order 
to maximize the knowledge gained about the underlying 
cosmological parameters from this high resolution experi- 
ment, it is important that our theoretical calculations are 
equally precise. Currently, uncertainty in the ionization 
history of the Universe remains one of the major factors 
limiting the accuracy of power spectrum calculations [e.g. 
see 0, y| . 

This fact has recently motived several studies on high 
precision computations of the cosmological hydrogen and 
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helium recombination history (see Sect. |TT] for a detailed 
overview) . The largest corrections thus far have been due 
to the non-equilibrium distribution of electrons within 
the hydrogen sub-states, the absorption of helium pho- 
tons by neutral hydrogen and the 2 3 Pi — l 1 So triplet- 
singlet transition in helium. Each of the additional phys- 
ical processes that have been discussed in the literature 
lead to > 0.1% level corrections of the ionization his- 
tory. While the individual changes partially cancel each 
other, the corresponding overall theoretical uncertainty 
in the CMB temperature and polarization power spectra, 
in particular at large I, still exceeds the level of 0.1%. It is 
also clear that some additional steps must still be taken, 
particularly in connection with the radiative transfer in 
hydrogen and the proper inclusion of two-photon pro- 
cesses (see Sect. Ill Cl for more details). These issues are 
still largely open questions and the subject of intensive 
ongoing study. 

However, the simultaneous and self-consistent inclu- 
sion of all these processes makes the computations very 
difficult and time-consuming. It involves the simulta- 
neous evolution of hundreds or thousands of differential 
equations for the occupation numbers of the individual 
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levels of the hydrogen and helium atoms. As an accurate 
calculation of the ionization history of a single cosmo- 
logical model using the current version of our full multi- 
level recombination code takes many hours or even days 
of computational time, standard parameter estimation 
methods require some approximation of the ionization 
history that is fast to evaluate. Currently the standard 
method used to evaluate the ionization fraction is the 
Recfast code Q, which models the hydrogen and he- 
lium atoms as effective 3-level systems. In order to ap- 
proximate the computations done using their multi-level 
recombination code Q, the authors introduced a fudge 
factor, /h, to artificially speed up hydrogen recombina- 
tion by about 14%. More recently, Wong et al. [f| up- 
dated Recfast to include a second fudge factor, &hc, to 
modify the recombination of helium. This change, which 
allows speeding up He I recombination, was motivated by 
the results of Switzer and Hirata Q, Kholupenko et al. 
[H and Rubiho-Martm et al. Q. It is intended to include 
the effect of the hydrogen continuum opacity during neu- 
tral helium recombination (see Sect. Ill B~2l) . 

Here we provide a different approach to computing a 
fast and accurate approximation of the recombination 
histor y. Our method uses a regression code based on 
Pico [13, IH, which we will henceforth refer to as RiCO, 
to model the ionization history, X e (z) = N c (z)/Nn(z), 
as a function of the cosmological parameters. Here N c 
denotes the number density of free electrons, and Nn 
is the number density of hydrogen nuclei in our Uni- 
verse. Since X e is a smooth function, the polynomial 
approximation used by RiCO is extremely accurate, and 
one retains nearly all of the information contained in the 
multi-level code used for training, the training code. 

As the accuracy of the training code is the major de- 
terminant of the accuracy of Rico's approximation, in 
Sect. [II] we review the physical processes included in our 
training code as well as an assessment of the processes 
that must still be considered and are currently under in- 
vestigation. In Sect. HTT1 we compare the ionization his- 
tory and corresponding power spectra computed with the 
multi-level recombination code used for the training of 
Rico to the different versions of Recfast. Also we in- 
clude a short discussion of how one could modify the ap- 
proach employed by Recfast to capture the new physics 
included in our full code. In Sect. IIVI we show that the 
simple regression scheme used by RiCO can accurately 
model the current full multi-level recombination calcula- 
tion, giving power spectra that are sufficiently accurate 
for even a cosmic variance limited experiment. 

Although the training code discussed here does take 
into account most of the important corrections discussed 
in the literature so far, it does not yet solve the prob- 
lem completely. Still more physical processes should be 
included and the results validated by independent codes. 
As RiCO can be trained equally well on any recombina- 
tion code it can facilitate crosschecks between those codes 
as well as a study of how approximations in these codes 
propagate to cosmological parameter constraints. 



Corrections to the Ionization History during Hydrogen Recombination 




FIG. 1: Corrections to the hydrogen recombination history. 

For the curves labeled '/-splitting' and 'our most recent 
100-shell computation' we compared with the RECFAST-code 
vl.3 after removing all of the switches (see Sect. IIII B[l , In 
the other cases we compared with the solution obtained with 
our code only including the non-equilibrium effects. 



Our target is to ensure that the ionization history can 
ultimately be calculated to sufficient precision to avoid 
biasing parameter estimation from the next generation of 
CMB experiments. By making it easy to propagate ad- 
vances in the calculation of the ionization history through 
to predictions of the CMB power spectra with RiCO, fu- 
ture development can focus on the physics of recombina- 
tion and to a lesser degree on the computational efficiency 
of the training code. 

RiCO and its future updates will be made available at 
http : / / cosmos . astro . uiuc . edu/rico. 



II. DESCRIPTION OF THE MULTI-LEVEL 
RECOMBINATION CODE USED FOR 
TRAINING 

In this section we discuss the previously neglected 
physical processes that are accounted for in the multi- 
level recombination code 1 which here is used to train 
RiCO. As wc will show, the addition of these processes 
lead to small but significant changes in the CMB power 
spectra (see Sect. IIIip . A comparison of the impact of 
all the included processes on the recombination history 
is illustrated in Fig. [1] and [2] In Sect. Ill CI we also give 



1 This code was developed as an extension of the works of Rubino- 
Martm et al. and Chluba et al. [H. 
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an overview and short discussion of processes that will 
be included in a future update of RiCO. 



A. Processes included during Hydrogen 
recombination 

1. N on- equilibrium populations of angular momentum 
sub-states 

While the previous work of Seager et al. [B[ evolved the 
individual energy levels of Hydrogen, thereby allowing 
departures from Saha-equilibrium, it was assumed that 
for principal quantum number n > 2 the angular momen- 
tum sub-states were in full statistical equilibrium with 
each other. Generally, equilibrium between I sub-states 
is maintained through collisional interactions. During 
recombination, however, collisions are much weaker than 
radiative processes, lea ding to a departure from full sta- 
tistical equilibrium [l2|, Il3| 

In the computations of Chluba et al. [H[ this lead to 
a ~ 5 — 10% increase in the free electron fraction, X e , 
at very low redshift [z < 400). Furthermore, a ~ 0.6% 
decrease in X e at decoupling was found in comparison 
with Recfast. As pointed out in their work this latter 
correction was not connected with the departures from 
full statistical equilibrium but has been interpreted as a 
hint towards limitations of the effective 3-level approach 
used in Recfast . 

For hydrogen here we follow the evolution of up to 
100 /-resolved shells. We also include the /-changing col- 
lisional rates, however their impact on the recombina- 
tion history is very small (< 0.01% with a maximum at 
z ~ 900). This problem involves the simultaneous evo- 
lution of up to ~ 5000 very stiff differential equations; a 
task that is numerically very cumbersome even on mod- 
ern machines. The detailed treatment of this process is 
one of the most difficult and time-consuming aspects of 
the whole training code. 

Here we have now also updated all the physical con- 
stants according to the NiST-database 2 . We then found 
that the 0.6% decrease in X e at decoupling practically 
disappeared. This is mainly because previously the ef- 
fect of the reduced mass of hydrogen, which leads to a 
~ 0.05% correction of the ionization potential, was ne- 
glected in the code of Chluba et al. ^3j . 

We have also improved the evaluation of the photoion- 
ization and recombination rates, which for many levels 
takes rather long. In Chluba et al. [HI the recombina- 
tion rates for all levels were tabulated before the actual 
computation and detailed balance was used to infer the 
photoionization rates each time the system is evaluated. 
This treatment is possible as long as the photon and 
electron temperatures do not depart significantly from 



each other, but becomes less accurate at low redshifts 
(z < 800). Here we use the more general procedure de- 
scribed in Chluba and Sunyaev [14j and found excellent 
agreement with the full computation, but at significantly 
lower computational cost. 

Furthermore, we also discovered a numerical inaccu- 
racy in connection with the computation of the recom- 
bination rates at low redshifts, present in the code of 
Chluba et al. [l3| . Fixing this problem decreases the final 
correction at z < 600 — 800 by a factor of ~ 2. However, 
this change is not very important in connection with the 
CMB power spectra, and also the results for the hydrogen 
recombination spectrum remain practically unaltered. 

The final result for the changes in the ionization frac- 
tion connected with the departures from full statistical 
equilibrium is shown in Fig. [1] (curve labeled '/-splitting'). 
To probe the robustness of this result, we carried out 
two sets of comparisons with the program developed by 
Rubino-Martm et al. [l2| • Already in the initial version 
of this code the most recent physical constants from the 
NiST-database and also the effect of the reduced mass 
were included. We explored the robustness of the low- 
redshift behavior by using detailed computations with 
n ~ 3 hydrogen shells down to redshifts z > 200. We 
found agreement at the level of < 0.001%. As a second 
set of tests, we focused on the redshift interval around the 
peak of the visibility function, and we compared compu- 
tations with n = 5, 10, 15 and 20 hydrogen shells. In 
all cases the agreement was better than 0.001% at the 
redshift range between z ~ 900 and z ~ 1600, where the 
analysis was carried out. However, for lower redshifts and 
larger values of n it will still be very important to have 
independent confirmation of the result presented here. 

One should also mention that below redshifts z ~ 200 
we do not follow the full system anymore, but extend 
the solution to z = using a simple 3-level atom ap- 
proach similar to Recfast. In order to smoothly con- 
nect the solution we re-scale the derivative of the electron 
fraction using the information from the last point of the 
full computation. This approximation should have little 
effect on CMB anisotropy power spectra. Also the ac- 
curate treatment of recombination in this redshift range 
should still include the details of primordial chemistry 
(see Sect. IE}- 



2. Induced decay of the 2s-level of Hydrogen 

The 2s — > Is two photon transition plays an important 
role in the recombination of hydrogen [TBI, [r| as it pro- 
vides one of the primary channels for creating neutral hy- 
drogen 3 . This decay rate is generally computed assuming 
no background photon field. However, the background of 



2 http://www.nist.gov/, May 2008. 



3 About 57% of all hydrogen atoms became neutral via this route, 
and only ~ 43% through the Lyman-ct line [i7ll 
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CMB photons in the Universe leads to stimulated transi- 
tions of the 2s-state to the ground level. Assuming that 
the CMB is given by a pure blackbody with temperature 
T 7 = 2.725 (I + 2) K, it was shown that the induced emis- 
sion leads to a ~ 1% increase in the 2s — > Is transition 
rate [HI . This increase has the effect of speeding up hy- 
drogen recombination leading to a maximum change in 
the free electron fraction of about 1.3% at z ~ 1050 (see 
Fig. [TJ curve labeled 'induced 2s-decay'). 

Since only the pure CMB blackbody contributes signif- 
icantly to the stimulated 2s —>■ Is transition, it is possible 
to tabulate the effective transition rate as a function of 
temperature before the computation. However, for the 
inverse process (see Sect. Ill A 3|) the Lyman-a spectral 
distortion is very important, so we also include it here 
for consistency. This is accomplished by first running 
a 5-shell computation of the hydrogen recombination to 
obtain a sufficient estimate of the Lyman-a distortion 
within the considered cosmology. The effective 2s — * Is 
transition rate is then tabulated at 5000 redshift points 
and log-interpolated. The integrals were evaluated with 
relative accuracy e ~ 10~ 8 . 



3. Feedback of Lyman-a photons on the effective Is — > 2s 
absorption rate of Hydrogen 

Due to the large number of photons produced in 
Lyman-a transitions during hydrogen recombination 
there is a huge excess of radiation over the background 
spectrum in the far Wien tail of the CMB. As noted 
by Kholupenko and Ivanchik [l9[ , after some redshifting, 
these super-Planckian photons are able to re-excite elec- 
trons from the ground-state to the 2s-level, and therefore 
increase the effective Is — > 2s absorption rate. This de- 
crease in the rate of hydrogen recombination (see Fig. [TJ 
curve labeled Ts-2s Lyman-a feedback') practically can- 
cels the effect of stimulated 2s-decays (see Sect. Ill A2j) . 
leading to a net ~ 0.6% increase in the free electron frac- 
tion at z ~ 980 (see Fig. [TJ curve labeled 'induced 2s- 
decay + ls-2s Lyman-a feedback'). While the low red- 
shift behavior of this last result differs slightly from the 
curve given in Kholupenko and Ivanchik [191 ] , it was also 
recently obtained by Hirata [20l ]. 

As explained in Sect. Ill A 21 in order to include the 
Lyman-a feedback we run a 5-shell computation of hy- 
drogen recombination to obtain a sufficient estimate of 
the Lyman-a distortion and then tabulate the Is — > 2s 
absorption rate. We checked that using more shells for 
the simple run does not affect the results significantly. 



4- Feedback within the Hi Lyman-series 

Due to redshifting, all of the Lyman-series photons 
emitted in the transition of electrons from levels with 
n > 2 have to pass through the next lower-lying Ly- 
man transition, leading to additional feedback correc- 
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FIG. 2: Corrections to helium recombination history. For 
the inlay plot we compared the ionization history with the 
standard RECFAST-code vl.3 after removing all of the 
switches (see Sect. IIIIB[l . 



tions like in the case of Lyman-a absorption in the 2s-ls 
two-photon continuum [211 ] . However, in this case the 
photons connected with Lyn, arc completely absorbed by 
the Ly(n — 1) resonance and eventually all Lyn photons 
are converted into Lyman-a or 2s-ls two-photon decay 
quanta. This process delays hydrogen recombination, 
leading to a maximal correction to the ionization history 
of AN C /N C ~ 0.22% at z ~ 1050 (see Fig. [TJ curve la- 
beled 'Lyman-series feedback'). This result was recently 
confirmed by Kholupenko et al. [22T ]. 

To include this feedback we save the spectral distor- 
tion due to emission of photons in the Lyman-series up 
to some given nf ee d (typically we use nf ee( j ~ 5) after 
each time step. We then include the distortion in the 
evaluation of the excitation rate of the next lower-lying 
Lyman-series transition, assuming that all photons en- 
tering the resonance are absorbed. Due to the huge op- 
tical depth in the Lyman-series this procedure is well- 
justified. In this way the Lyman-series feedback always 
works like Lyn — > Ly(n — 1). Note that due to the dis- 
tance between neighboring Lyman-series resonances, the 
feedback occurs after the time it takes to redshift into 
the next transition. 



B. Processes included during Helium 
recombination 

For He n-recombination a list of processes similar to 
hydrogen can be formulated. In particular C. Hirata and 
E. Switzer have rigorously studied these and several ad- 
ditional processes 0, [H, 24]. However, since the details 
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of helium recombination history are not strongly propa- 
gating to the computation of the CMB power spectra, as 
a first step here wc only include the two dominant cor- 
rections, which actually are due to processes that have 
no analog in the case of hydrogen recombination. Ad- 
ditional corrections (see discussion in Sect. Ill Cj) will be 
taken into account for a later release of RlCO. 

One should also mention that for the computation in- 
volving helium wc only include 5 j-resolved shells for the 
singlet and triplet states [for more detail see [§]. We 
found that within the current numerical approach and 
related approximations this leads to sufficient precision 
(see Sect. lITCTj) . 

1. Spin-forbidden Hei-2 3 Pi — l 1 So transition 

Like for the Lyman-a photons in hydrogen, due to the 
low expansion rate of the Universe, He i-2 1 Pi — l 1 So pho- 
tons have a very hard time escaping from the interac- 
tion with neutral helium during He ii — > He i recombi- 
nation. Therefore also the recombination of helium is 
strongly delayed as compared to the Saha case, and the 
Hei-2 1 S - 1 1 S two-photon transition [A 2l ~ 51.4s -1 ; 
[25j plays a very important role in defining the rate of 
He i-rccombination. Within the standard computation 
the Hei-2 1 So — l 1 So two-photon decay channel allowed 
~ 31% of helium atoms to recombine, while the Hei- 
2 1 Pi — l 1 So transition contributed ~ 69% [e.g. see table 

in Eg. 

Dubrovich and Grachev E3] realized that in addition 
to the He 1-2 1 So — l 1 So two-photon transition, the spin- 
forbidden Hei-2 3 Pi — l 1 So transition helps helium to re- 
combine even th oug h it also has a very low probability 
[A ~ 1 77 s _1 ;EE E3 as compared to the He 1-2^1 -1 1 S 
resonance (A ~ 1.80 x 10 s ). This process, which 
was not included in the original version of Recfast Q , 
leads to a ~ 1.1% decrease in the free electron fraction 
at z ~ 1750. Note that in terms of the helium recombi- 
nation history this corresponds to a very large effect (see 
Fig.©. 

By including the spin-forbidden He i-2 3 Pi — l 1 So tran- 
sitions into the recombination code one finds that about 
~ 17% of all helium atom became neutral through the 
Hei-2 1 So — l 1 So two-photon transition, ~ 40% via the 
He i-2 1 Pi - 1 1 S transition and ~ 43% through this 
new channel [26| . The comparison shows that the spin- 
forbidden Hei-2 3 Pi — l 1 So transition is actually one of 
the dominant channels for helium recombination. This 
conclusion is true even when the hydrogen continuum 
opacity-effect (see section Hi B 2[) is also included. In this 
case only 8% of all helium atoms recombine via the He I- 
2 1 So — l 1 So two-photon decay, while the rest is shared 
between the Hei-2 1 Pi — l 1 So and spin-forbidden Hei- 
2 3 Pi — 1 So transition, with a ~ 2% contribution from 
other direct transitions to the ground-state 0. 

This process can be self-consistcntly included by sim- 
ply adding another term similar to the normal resonance 



transitions to the set of rate equations [27j, a modifica- 
tion that is now also accounted for in Recfast vl.4 [||. 
Note that for this transition the Sobolev-escape proba- 
bility should be taken into account, since the Sobolev 
optical depth still reaches ~ 3 close to its maximum. 



2. Absorption of He I photons by neutral hydrogen 

In addition, the recombination of neutral helium is 
sped up due to the absorption of Hei-2 1 Pi — l 1 So and 
2 3 Pi — l 1 So-photons by the tiny fraction of neutral hy- 
drogen that is already present at redshifts z < 2400. 
While photons arc emitted close to optically thick Hei- 
resonanccs, a small part of them can be absorbed in the 
Lyman-continuum of hydrogen, allowing additional elec- 
trons to settle to the ground state of helium. This pro- 
cess was mentioned by P. J. E. Peebles in the mid 90's 
[see remark in Q , but only recently taken into account 
by Switzer and Hirata [23| and others [1, Q . It leads to 
the largest correction during the epoch of Heii — > Hei- 
recombination that has been investigated so far in the 
literature, strongly accelerating the recombination of he- 
lium below z ~ 2000 (sec Fig. [2 thick solid line). 

Since this problem is connected with important de- 
tails in the radiative transfer of optically thick Hei- 
2 : Pi — l : So and 2 3 Pi — l 1 So, it is a very time-consu ming 
task to solve self-consistently. Switzer and Hirata [23[ 
studied this problem using a Monte-Carlo approach to 
solve the quasi-stationary line-transfer problem. They 
tabulated the escape probability for different values of 
the Sobolev-optical depth and hydrogen continuum opac- 
ity and then interpolated these in the actual recombina- 
tion calculation. Independently, Kholupenko et al. [1] 
investigated this problem using a simplified analytical 
approach, which slightly underestimated the acceleration 
of helium recombination. They provided simple expres- 
sions for the correction to the Sobolev escape probability 
of the He i-2 1 Pi - 1% and 2 3 Pi - 1%, which were 
later also used by Wong ct al. @ to fudge the helium re- 
combination history. Rubiho-Martm et al. @ confirmed 
the results of Switzer and Hirata [23[ using the results 
of a Fokker-Planck approach to solve the line-transfer 
problem [jjcj. They also gave a simple ID-integral which 
neglects the redistribution of photons by resonance scat- 
tering but still reproduces the correction to the escape 
probability rather well. Very recently Kholupenko et al. 
[3lj reconsidered this problem in more detail analytically 
and found good agreement with Switzer and Hirata [231 ] 
and Rubiho-Martm et al. Q. 

In this work we follow the approach of Rubiho-Martm 
et al. Q and include the acceleration in the helium recom- 
bination history using their ID-integral expression [see 
Eq. (B.3) in©] for the correction to the escape probabil- 
ity. Here we do not include the redshift dependent fudge- 
function that was introduced by Rubino-Martm et al. [9j 
to account for the additional acceleration because of par- 
tial redistribution. In the future we plan to solve the full 
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time-dependent problem including the radiative transfer 
in the optically thick lines, a task that will also become 
necessary for hydrogen recombination. 

C. Additional processes that will be considered in 
the future 

There arc a number of additional processes, both for 
hydrogen and helium recombination, that have already 
been addressed in the literature, but were not taken 
into account in the training code. As mentioned above, 
C. Hirata and E. Switzer have rigorously studied several 
additional subtle processes (e.g. feedback; two-photon 
processes; effect of electron scattering; isotope shift of 
3 He) within the context of helium recombination [see 
0, HE H- Also for hydro gen recombination additional 
processes have been discussed, however not all of these 
studies are concluded yet. Some of these additional pro- 
cesses lead to rather small additional corrections, but 
others may still be important. Below we give a short 
overview about some of the work that must still be done. 
We are planning to make a careful survey in the near 
future. 



1. Incompleteness of the atomic models 

At low redshifts (z < 800) the rate of hydrogen recom- 
bination is strongly controlled by the effective recombina- 
tion coefficient, which itself depends on the completeness 
of the atomic model of hydrogen. Here we included only 
up to 100 shells in our computations, but as pointed out 
earlier (l3j . in terms of the ionization history this may 
still not be enough, and extension to ~ 200 — 300 shells 
may be required. However, for the computation of the 
CMB power spectra these corrections probably are not 
very important, as even with 75 shells rather converged 
results seem to be obtained (see Sect. IIV[) . Pushing to a 
larger number of ^-resolved shells is not trivial, and also 
more accurate collisional rates may be required. 

As mentioned above, for the computation involving he- 
lium we only include 5 j -resolved shells. This approxima- 
tion seems to be sufficient at the current level of precision 
since, unlike the case of hydrogen recombination, the dy- 
namics of helium recombination are much less controlled 
by the effective recombination rate, which is strongly con- 
nected with the completeness of the atomic model and 
can be fudged to some level. For helium recombination 
the escape of photons is much more crucial. 

2. Escape of H I Lyman-ct photons 

The escape of photons from the optically thick H I 
Lyman-a line is usually modeled using the Sobolev- 
approximation. The validity of this approximation dur- 
ing hydrogen recombination has been investigated several 



times [e.g.|H, [Hjj, but at the one percent-level a full con- 
firmation is still necessary. For example, Grachev and 
Dubrovich [34| recently claimed that the line-recoil effect 
leads to a ~ 1.3% speed-up of hydrogen recombination. 
Within the quasi-stationary approach this results seems 
robust ■ Given the importance of this problem to the 
calculation of the power-spectra, independent checks are 
necessary. This is now being investigated in detail [sod - 



s'. Two-photon transitions from higher levels 

Along with the 2s — > Is two-photon transition there 
are also allowed two-photon transitions from higher levels 
to the ground state. Within the context of hydrogen and 
helium recombination these transitions were first studied 
by Dubrovich and Grachev [27J . They predicted a ~ 
5% decrease in the free electron fraction at z ~ 1200, 
however in the computations of the effective two-photon 
decay rates of the ns and nd-levels they only included 
the first non-resonant term (i.e. due to the dipole matrix 
element connecting ns/nd — ► np) in the infinite sum over 
intermediate states. 

Using rate coefficients for the two-photon decay of the 
3s and 3d-levels in hydrogen as computed by Cresser 
et al. [35j . Wong and Scott [26[ found that Dubrovich 
and Grachev [27| overestimated the impact of two-photon 
transitions on the ionization history by about one order 
of magnitude. However, the calculation of Cresser et al. 
[35l ] was incomplete. For example they did not include 
the largest non-resonant term (due to the dipole matrix 
element connecting ns/nd — > np) in their calculations 
[3^ . Also physically it is very difficult, if not impossible, 
to separate the 'pure' two-photon decay rate from the 
resonant contributions [see discussions in [2(| [H, H3| , e.g. 
because of non-classical interference effects. 

Later this problem was reinvestigated in more detail, 
and a lower limit for the impact of two-photon decays 
during hydrogen recombination was derived, implying 
that a decrease of slightly more than 0.3 — 0.5% in the 
free electron fraction at z ~ 1150 can still be expected 
[36l | . This estimate was obtained by taking into account 
departures of the full two-photon line profiles from the 
Lorentzian shape in the very distant, optically thin part 
of the red wing of the Lyman-a line. According to these 
computations, the two-photon decays from s-states seem 
to slow hydrogen recombination down, while those from 
d-states speed it up. In addition it was shown that the 
slight net acceleration of hydrogen recombination seems 
to be dominated by the 3s and 3d contribution (36j . 

However, it was pointed out that the final answer can 
be only given using a full radiative transfer computation, 
which also takes into account the effects of partial fre- 
quency redistribution and the feedback of photons from 
the blue side of the Lyman-ct resonance [e.g. see |36|. 
Very recently Hirata [2(| showed that including these as- 
pects of the problem, along with the induced 2s-decay 
and the feedback of Lyman-a photons on the effective 
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Is — > 2s absorption rate, modifications in the ionization 
history of the order of ±1.3% can be expected. Also 
2s — > Is Raman-scattering seems to play some role. This 
would be a very important effect and we are currently 
investigating it. 

Also during helium recombination two-photon transi- 
tions and Raman-scattering are important. Hirata and 
Switzer [24| investigated the effect of these, as well as the 
stimulated 2 1 So — > l 1 So two-photon decay, feedback by 
photons from the Hei-2 1 Pi — l x So and spin-forbidden 
Hei-2 3 Pi — l 1 So transition on the l^So — > 2 1 So absorp- 
tion rate. However, the overall changes due to the two- 
photon transitions and Raman-scattering seem to be in- 
significant [24j]. 
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4- Feedback due to helium lines 

The feedback of high frequency photons released dur- 
ing helium recombination should also affect the dynam- 
ics of hydrogen recombination. Here it is interesting that 
most of the high energy photons from He in — > Heii 
will be reprocessed by neutral helium before they can 
directly affect hydrogen. Since the number of high fre- 
quency 7(Heii) photons is comparable with the number 
of helium atoms this should still have a rather strong ef- 
fect. Note that although the ionization history during 
He in — > He ii-recombination has a very small impact 
on the CMB power-spectra, the exact time-dependence 
of the 7(Heii) — > Hei feedback is related to how fast 
He in — > He ii-recombination occurred, so that it may 
still deserve careful investigation. 

Also feedback of the Hei line during helium recombi- 
nation still has to be included, and was shown to have 
a notable delaying effect [23[ . However in particular the 
feedback of Hei-2 1 Pi - 1% photons on the 2 3 Pi - 1% 
transition may in addition require a fully time-dependent 
treatment, since both resonances are only separated by 
~ 1% in frequency, or ~ 600 Dopplcr widths relative to 
the Hei-2 1 Pi - l^So resonance at z ~ 2500. 

Finally, those photons emitted by neutral helium can 
directly feedback on hydrogen, but in order to take this 
feedback into account more detailed computations of the 
helium recombination spectrum are required. This pro- 
cess could still affect the hydrogen recombination history 
on a level exceeding 0.1%. 



5. Details of the primordial chemistry at low redshifts 

For the initial computations in connection with Rec- 
fast [j| some aspects of the primordial chemistry were 
included. These have not yet been taken into account in 
our code, but are expected to have only a small effect, 
mainly at redshifts z < 200. However, it is in principle 
easy to include these, as well as extending the chemical 
network using updated rate coefficients [e.g. see[H|. 



FIG. 3: Comparison of the angular power spectra computed 
using the ionization history of our most recent 100 shell code 
versus the default settings in Recfast vl.3 and vl.4. The 
two black lines denote our performance benchmark of ±3/£ 
which is roughly 2 times cosmic variance for the special case 
of fitting a single amplitude parameter (see also Seljak et al. 

1). 



6. Other small corrections 

Similar to the effect of 3 He on helium recombination 
[7j , also deuterium should affect hydrogen recombination. 
Since due to the isotope shift the deuterium Lyman-a res- 
onance is on the blue side of the hydrogen Lyman-a line, 
this should slow down hydrogen recombination. How- 
ever, since the abundance of deuterium is so small, and 
because at z ~ 1100 the shift is only ~ 12 Doppler- width 
of the Lyman-a line, this process is probably negligible. 

For the feedback of the CMB spectral distortion gen- 
erated during recombination on the Is — > 2s-absorption 
rate we only took into account the hydrogen Lyman-a 
distortion. However, the high frequency part of the 2s- 
distortion itself may also have to be included, probably 
leading to another small correction to the ionization his- 
tory. 

III. COMPARISON WITH RECFAST 

The net effect of inclusion of the physics listed in 
Sects. Ill Al and IIIBI leads to an increase in the free elec- 
tron fraction during hydrogen recombination. Close to 
the maximum of the Thompson visibility function, which 
is most relevant for computations of the CMB power 
spectra, one can observe a 0.7% correction. At low z 
the error increases to ~ 3% at z ~ 200. During he- 
lium recombination the correction reaches about —3% at 
z ~ 1800 (see Fig. for details). 

In Figure [3] we show the corresponding change in the 
temperature and polarization power spectra when com- 
paring with the output of Recfast vl.3 and vl.4 us- 
ing the default settings. We used the publicly avail- 
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able code 4 CAMB 39j to compute the CMB power spec- 
tra. The cosmological model shown has parameter val- 
ues n B = 0.0444, fi M = 0.2678, Q K = 0, H = 71 
km/s/Mpc, T cmb = 2.725 K, and Y p = 0.24. The solid 
and dotted lines show the fractional error between our 
100 shell calculation and Recfast for temperature and 
polarization spectra respectively. 

The figure shows that the corrections during helium re- 
combination contribute about 1/2 to the average change 
in the CMB power spectra. Inclusion of the spin- 
forbidden Hei-2 3 Pi — l 1 So transitions and fudging of the 
Hei-2 1 Pi — 1 So escape probability to account for the ef- 
fect of the continuum opacity from neutral hydrogen as 
done in Recfast vl.4 leads to a factor of 2 improvement 
in the accuracy of the power spectra. The remaining 
error, which is due to inaccuracy in the model of hy- 
drogen used by Recfast, remains above the 0.5% level 
at I > 1000. This may not be sufficient for analysis of 
Planck data. 

Also plotted are dotted lines corresponding to ±3/£, 
as suggested by Seljak et al. [|| for a benchmark. These 
lines correspond to roughly 2 times cosmic variance for 
the special case of exploring the constraints on a sin- 
gle cosmological parameter. Since correlated errors over 
several I values can lead to bias in parameter constraints 
even if the error at individual £'s are less than cosmic 
variance, these lines effectively denote an estimate of the 
minimum error to which any experiment may be sensi- 
tive. The upcoming Planck satellite mission should probe 
the temperature power spectrum beyond I ~ 2500 and 
be cosmic variance limited out to I ~ 1500. Thus the ac- 
cumulated error in the power spectrum due to even sub- 
percent level errors in the ionization history may have a 
significant impact on parameter constraints derived from 
Planck data, in particular when considering the potential 
bias on inflation parameters (e.g. the spectral index and 
its running). 



A. Modification to the Recfast fudge factors 
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FIG. 4: Comparison of the ionization fraction during 
hydrogen (top) and helium (bottom) recombination between 
Recfast vl.4 and our training code. The ionization history 

is shown using the default fudge factors (/h = 1.14, 
fee = 0.86) and the value that minimizes the fractional error 
in X e (fn° — 1.106, 6hc = 0.85). The top panel also shows 
the value of /h that minimizes the cosmic variance weighted 
error in the power spectrum (f^ e = 1.065). We have also 
noted the artifacts in the ionization fraction due to the 
switches in Recfast (sec Sect. IIIIB|) . 



Some of the improvements to the calculation of the 
ionization history discussed in Sects. Ill Al and IIIBI can 
be captured in Recfast by modifying the hydrogen and 
helium fudge factors, denoted by /h and 6hc respectively. 
For hydrogen, adjusting /h to minimize the fractional er- 
ror in the ionization history against our most recent 100 
shell calculation over the redshift range 600 < z < 1200 
gives an optimal value of = 1.106. This number is 
comparable to the one obtained in Rubino-Martm et al. 
[l2l |. where they found = 1.10 for 30-shells, and is 
roughly a 3% decrease compared to the default value of 
/h = 1.14 used in Recfast. We have performed the 
same optimization over 20 other models and find the op- 



4 http://camb.info 



timal value of /h to lie in the interval [1.102,1.107], in- 
dicating that there is a only a small dependence of the 
optimal fudge factor on the cosmological parameters. We 
note here that, as can be inferred in Fig[TJ the derived 
value for /h is particularly sensitive to the lower limit of 
the considered redshift interval. 

A similar optimization of the helium fudge factor over 
redshifts 1400 < z < 3200 gives 6 Ho = 0.85, a small 
change compared to the default value of 0.86 used in 
Recfast vl.4. The modification to the ionization his- 
tory from Recfast by adjusting these two fudge factors 
is shown in Figure [5] One should also mention that some 
differences are expected since the effect of partial redistri- 
bution was neglected here. This should not have a large 
impact but will be considered in a future release of RiCO. 

If instead the goal is to minimize the error in the power 
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FIG. 5: Fractional error in the temperature and polarization 

angular power spectra between Recfast vl.4 and our 
training code. The hydrogen fudge factor in Recfast is set 
as /h' = 1.065 which minimizes the cosmic variance 
weighted error in the power spectra. The two solid black 
lines denote our performance benchmark of ACe/Ce — 3/£ 
(see Seljak et al. 0]). 



spectra we find the optimal value of the hydrogen fudge 
factor to be / H * = 1.065. As shown in Figure[5j using the 
modified fudge factor captures some of the changes to the 
ionization history. In particular, using fg* as the fudge 
factor improves the fit in the region where the visibility 
function is large by sacrificing accuracy at z < 800. This 
gives a factor of 2 improvement in the power spectra (see 
Figure [5]) . The further change in helium recombination 
between our code and Recfast vl.4 has only a tiny effect 
on the power spectra. 



B. Effect of the switches in Recfast 

In order to ease the solution of the effective 3-level sys- 
tem used in the Recfast code, some switches were intro- 
duced 0], that are still present in the current version 1.4. 
Starting at redshift ~ 2800 for the standard cosmology 
the RECFAST-code only solves the differential equation 
for helium, using the Saha solution for hydrogen. Then 
when the ratio N p /N-r < 0.99 (for the standard cosmol- 
ogy this occurs at z ~ 1550) also the hydrogen equation 
is solved. Depending on the derivative for the matter 
temperature, finally also the temperature equation is in- 
cluded to the full system (for the standard cosmology this 
occurs at z ~ 820). 

In Fig. [S] we show the effect of these switches, which we 
avoided in our version 5 of Recfast. For this we changed 



old values for G & mHe/mH 
new values for G & inHc/mH 
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FIG. 6: The effect of switches in Recfast vl.3 and vl.4. For 
Recfast vl.4 we have not included any of the recent helium 
corrections. In both cases we compared with our version of 
Recfast, which avoids the switches (see Sect. HUB)) . At 
high redshifts one can see the small effect of mn e /irm. 



to a solver for stiff differential equations from the Nag- 
Library 6 . The most relevant effect is the switch for the 
temperature equation, which is producing a deviation of 
the order of ~ 0.2% at redshifts z ~ 900, with respect to 
our solution in which the temperature equation is solved 
at all redshifts. Note that in order to accurately sample 
the peaks that result from the switches as well as to main- 
tain consistency with CAMB we increased the sampling 
of redshift points used by Recfast to 10 4 . 

For completeness, Fig. [6] also presents the comparison 
with the case in which the corrections due to the new val- 
ues of G and mHe/wiHj pointed out by Wong and Scott 
[26| . are not included. These two modifications were in- 
troduced in Recfast vl.3. One can see that at z < 2000 
these lead to corrections that are below ~ 0.05%, and ac- 
tually only reach ~ 0.1% at z > 6000. This modification 
was already taken into account in the initial version of 
the work by Rubiho-Martm et al. [l2l |. 

Note that for the comparisons in this section, our 
"Recfast computation" does not take into account any 
of the helium corrections, which are incorporated to Rec- 
fast vl.4 using fudging. 



5 Note that also the small bug in the differential equation for hy- 
drogen present in Recfast vl.2 was corrected. This bug leads 



to a ~ 0.1% decrease in X a at z ~ 1450. 6 See |http: / /www. nag. co.uk/numeric/ 
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Ionization History in Comparison with Recfast 
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FIG. 7: Comparison of the ionization history of our most 
recent 100-shell code with Recfast vl.3 after removing all 

the switches (see Sect. IIH"B)) . Also included are the 
variations in Recfast vl.3&vl.4 for different values of the 
hydrogen fudge factor. 
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FIG. 8: Relative change in the number density of free 
electrons when allowing for a constant relative change in the 
Sobolev-escape probability of the H I Lyman-a transition as 
labeled. All curves where computed using Recfast. At 
maximum AN C /N C ~ -jAP/P. 



C. To what extent does fudging work? 

In this Section we discuss our perspective of how far 
one can go with fudging, and what kind of modification 
may help to further improve the result with a simplified 



approach similar to Recfast. However, to us 'fudging' 
at some level either becomes equivalent with fitting the 
whole problem, or the performance of the corresponding 
algorithm would probably also be insufficient for cosmo- 
logical parameter estimations. This is the reason why 
we decide to take a direct fitting approach (see Sect. IIV|) 
leading to Rico, instead of a detailed investigation of the 
possibilities for more sophisticated fudging. 



1. Fudging hydrogen recombination 

Physically the fudge factor /h introduced to speed up 
hydrogen recombination increases the effective recombi- 
nation coefficient to the second shell. This implies that 
the hydrogen fudge factor can only alter the recombina- 
tion history significantly when the recombination dynam- 
ics truly depend on this coefficient. However, as Fig. [7] 
illustrates this is only the case at low rcdshifts. While 
to some extent the differences with our full recombina- 
tion calculation can be reduced to optimizing the value 
of /h (see Sect. IIII A[) . clearly not all the differences can 
be erased (see Fig. 2] and [5]). In addition there is some 
slight dependence of /h on cosmology (see Sect. IIII Al) . 
which one may still be interested in. 

As proposed by Chluba et al. [l3j |. the next level of 
fudging is using a redshift dependent function for /h- 
Looking at Fig. [7] one can realize that in particular 
at z ~ 1100 this fudge-function should depend rather 
strongly on redshift, and that for the high redshift part 
this approach will likely not lead to very good results. 
Additionally the cosmology dependence of this fudge- 
function will probably be significant. In the end it is 
equivalent to directly fitting the correction for the stan- 
dard concordance model only. 

In order to capture some of the corrections from a 
more physical point of view, one could think about sim- 
ple extensions of the effective 3-level atom. The induced 
two-photon decay (see Sect. Ill A 2[) can be incorporated 
rather easily, since only the effective decay rate due to 
blackbody CMB photons should be replaced [ll[ . For 
the feedback of the Lyman-a distortion on the effective 
ls-2s absorption rate (see Sect . Ill A3)) one should provide 
an approximation for the Lyman-a distortion, which in 
principle can be done analytically [l^ |. Including these 
two corrections should therefore be possible. 

In order to include the Lyman-series feedback (see 
Sect. Ill A~4|) one could simply take more shells, say 5 — 10, 
into account, and hence model this process more or less 
sclf-consistently. However, runs with 5 shells will take 
longer, and one should provide effective recombination 
coefficients for all the additional shells in order to cor- 
rectly model the Lyman-series lines. An alternative an- 
alytic treatment was recently presented by Kholupenko 
et al. [p]. 

In order to include changes close to decoupling one 
could also modify the H I Lyman-a escape rate, or alter- 
natively the effective two-photon decay rate of the the 2s 
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level. As an example, in Fig. [8] we illustrate how changes 
in the H I Lyman-a escape affect the ionization history. 
This could certainly help to take into account some more 
details of the radiative transfer problem, which are still 
under discussion and not included here (see Sect. Ill C2|) . 
Also the effect of two-photon processes from high lev- 
els (see Sect. Ill C 3|) and to some extent the details of 
radiative transfer (see Sect. Ill C~2|) could probably be ap- 
proximated within such an approach. However, this will 
certainly require a carefully calibrated redshift dependent 
fudge-function, as was also used by Rubiho-Martin et al. 
Q in the case of helium recombination. Again one then 
should check the cosmology dependence of this fudge- 
function, which, depending on the desired accuracy, may 
vary significantly. It is clear that within such an ap- 
proach the connection to the real physical processes will 
no longer be obvious. 



2. Fudging helium recombination 

For helium recombination introducing a fudge-factor 
for the effective helium recombination coefficient would 
not help at all. Even if one increases the recombina- 
tion coefficient by more than a factor of two, the he- 
lium recombination history basically docs not change. 
Therefore one directly has to fudge the escape rate of 
Hei-2 1 Pi - 1% and 2 3 Pi - 1% photons, as it has 
been done recently by Wong et al. [6[. Again, depend- 
ing on the desired accuracy, one should allow for a red- 
shift dependent fudge- function as was already proposed 
by Rubiho-Martin et al. @, since only the low redshift 
tail of helium recombination is accurately reproduced by 
Recfast vl.4. However, since details in the helium re- 
combination history are not strongly propagating to the 
computation of the temperature and polarization power 
spectra, fudging is probably sufficient from this point of 
view. Still, using Rico would enable one to reproduce 
the helium computations at a much higher level of accu- 
racy and with significantly less effort, even when more 
detailed computations become available. 



IV. APPROXIMATING THE FULL 
RECOMBINATION CALCULATION 



Instead of introducing new fudge factors or modify- 
ing existing ones to reproduce the results of the full re- 
combination code using the simplified 3-level model, we 
propose directly fitting to the results of the accurate cal- 
culation using a training set of ionization histories with 
a regression code similar to Pico [HJIIIj]- Since the ion- 
ization history is a smooth function of the cosmological 
parameters, RiCO is able to accurately capture most of 
the variation in this function. 

Ignoring massive neutrinos, the ionization history de- 
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FIG. 9: Comparison of the ionization history computed by 
our code for varying numbers of shells. The fractional error 
is taken against the 120 shell calculation. Note that little 
physics is lost by using the 100 shell approximation 
compared to the 120 shell and that the 60 shell computation 
remains below our minimum cosmic variance estimate until 

around I ~ 2000. The two solid black lines denote our 
performance benchmark of ACe/Ct = 3/1 (see also Seljak 
et al. 



pends on 6 cosmological parameters 7 : the baryon density 
fifth 2 , the dark matter density ^cdm^ 2 , the curvature 
density fix, the Hubble constant Hq, the temperature of 
the CMB T cm b and the helium mass fraction Y p . The 
training set for RiCO is generated based on constraints 
from the WMAP 3 year results It is therefore con- 
venient to use 9, the ratio of the sound horizon to the 
angular diameter distance at decoupling, instead of H as 
there will be significantly less correlation in the parame- 
ters [!l[. The parameters ^b^ 2 , Hcdm^ 2 , and 9 are 



7 If necessary one could probably reduce the number of indepen- 
dent variables by changing to a parameterization that encodes 
the variables that directly impact the recombination calculation. 
However currently this is not a limiting factor. 



12 




le _ 06 U , 1 , , , 1 1 1 1 1 , , , 1 , , 1 1 1 1 , 1 u 

500 1000 1500 2000 2500 3000 

z 

FIG. 10: Performance of RlCO in computing the ionization 
history based on 500 training cases computed with 60 shells. 
This test was done over 400 separate test cases. The lines 
denote the mean fractional error along with the fractional 
error that bounds 99% of the test cases. 



chosen to lie within 25 log-likelihoods of the peak of the 
WMAP 3 year likelihood. A method for efficiently find- 
ing these points is discussed in Fendt and Wandelt [111 ]. 
Values for the current temperature of the CMB and the 
helium fraction are chosen uniformly from the intervals 
[2.72,2.73] H3] and [0.232,0.258] jH respectively. In the 
parameter estimation framework, outside of this region 
we just use the standard Recfast code. However, these 
points will have low likelihood and should not effect the 
results the results of parameter analysis. 

We train RlCO using 500 parameter sets and corre- 
sponding ionization histories. For this work, we chose to 
use 60 shell models as they remain very accurate out to 
£ ~ 2000. This is shown in Figure [9] where we have plot- 
ted the fractional error using various numbers of shells 
compared to the 120 model. The two black lines again 
correspond to ±3/£. Note the small difference between 
the 100 and 120 shell model, indicating that there is little 
additional information gained about the power spectra 
from tracking shell populations beyond 100. Also, the 60 
shell line remains below our strict error tolerance until 
roughly £ ~ 2000 indicating that Planck will not be sen- 
sitive to the error from using the lower shell calculations. 

To test the accuracy of Rico we generated a test set 
of 400 models chosen in the same manner as the training 
set. The models used for this test are not included in 
the training set. As in the training set, the recombina- 
tion histories that make up the test set were run using 
60 shells in order to access only the error incurred by 
approximating the ionization history with RlCO. Figure 
1101 shows the error in the ionization history computed 
by RlCO against the full recombination code. At high z, 
there is little variation in the training set so RlCO has 
no trouble approximating the ionization fraction. After 



z ~ 1700, as He I recombination begins, until the end 
of hydrogen recombination the error from RlCO remains 
below 0.03% for 99% of the test cases. The correspond- 
ing error in the power spectra from using RlCO is thus 
negligible (see Fig. [TTj) . 

This demonstrates that RlCO essentially propagates all 
of the information from the ionization history into the 
power spectra. The ability to accurately compute the 
Ce's is primarily limited by the accuracy of the full recom- 
bination code used to train RlCO. Since RlCO uses the 
same fitting methodology as Pico it is possible to extend 
the parameter space or add additional model parameters 
(for example parameterizing uncertainties in the recom- 
bination calculation). This would increase the one-time 
training cost, but would not lead to an appreciable addi- 
tional cost in the evaluation of the ionization history. If 
necessary, the accuracy of RlCO can be further increased 
by adding more points to the training set as well as by 
using a higher order regression. Again, the only compu- 
tational penalty is to the one-time training step. 

While the downside of this method is the requirement 
of running the full code ~ 500 times to generate a train- 
ing set, this can be done completely in parallel. Large 
computing clusters or distributed computing projects are 
ideal for exactly this type of application. Also, this train- 
ing cost does not affect the user of the code. The advan- 
tage is that there is no need to find optimal approxima- 
tions using a simple physical system. RlCO can just be 
trained using results from the most accurate code avail- 
able. Also, RlCO is trained to compute the ionization 
history over a volume of parameter space and not sim- 
ply optimized based on a single model at the peak of the 
parameter likelihood. 



V. CONCLUSION 

We have presented a new code designed to compute 
the ionization history of the Universe. The code includes 
previously neglected physics that leads to changes in the 
ionization fraction at the 2 — 3% level in some redshift 
regions. This change leads to a correction to the CMB 
power spectra of more than 1% at £ ~ 3000. As uncer- 
tainty in the ionization fraction is the main contributor 
to error in the theoretical CMB anisotropy power spec- 
tra, this code represents a significant step in our ability 
to compute the power spectra to a significant precision 
for upcoming experiments. 

While it is possible to attempt to capture the changes 
to the ionization fraction by modification of the hydro- 
gen fudge factor in Recfast, some residual error remains 
near the peak of the visibility function. The problem of 
correctly introducing fudge factors can be avoided en- 
tirely by using a regression code based on a training set 
of cosmological parameters and their corresponding ion- 
ization histories. A simple polynomial fit is sufficient to 
compute ionization fraction at the level of ~ 0.01% over 
the volume of parameter space relevant to current exper- 
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FIG. 11: Fractional error in the power spectra using the 
recombination history computed by RlCO. The plot includes 
200 models from our test set and the comparison has been 

made with the power spectra using the recombination 
history from the 60 shell run of our full recombination code. 
The two solid black lines denote our performance benchmark 
of ACe/d = 3/1 (see Seljak et al. 0). 



imental data. 

As it is straightforward to train RlCO on multiple re- 
combination codes, our algorithm enables the propaga- 



tion of approximations made to the ionization history 
to biases on the cosmological parameters when analyz- 
ing data. Since implementing the many new physical 
processes already discussed in the literature in an ac- 
curate and robust manner is very challenging, we hope 
that this ability facilitates a model by model cross val- 
idation among the recombination codes that have been 
developed. Given the complexity of the physical pro- 
cesses involved in the computation of the cosmological 
recombination history, such a comparison will become 
very important in order to ensure that the final result 
is correct. Our goal is to ensure that the ionization his- 
tory can ultimately be calculated to sufficient precision to 
avoid biasing parameter estimation from the next gener- 
ation of CMB experiments. By making it easy to propa- 
gate advances in the calculation of the ionization history 
through to predictions of the CMB power spectra with 
RlCO, future development can focus on the physics of re- 
combination and to a lesser degree on the computational 
efficiency of the physical recombination code. 

RlCO and its future updates will be made available at 
http : //cosmos . astro . uiuc . edu/rico. 
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